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HIGHLIGHTS 


►  The  first  detailed  model  of  a  direct-borohydride/oxygen  fuel  cell  is  developed. 

►  The  performance  is  simulated  and  compared  to  experimental  results. 

►  The  performance  on  Pt  and  Ni  anodes  is  simulated  and  the  differences  are  explained. 

►  Extensions  to  the  model  are  outlined. 

►  Key  property  values  missing  from  the  open  literature  are  highlighted. 
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A  detailed  mathematical  model  for  the  direct  borohydride/02  fuel  cell  is  developed.  The  activation 
polarizations,  mass-transport  limitations  and  resistances  to  charge  transport  in  each  component  of  the 
cell  are  explicitly  incorporated.  The  anode  kinetic  mechanism  is  based  on  direct  borohydride  oxidation, 
borohydride  hydrolysis  and  the  full  Tafel-Volmer-Heyrovsky  mechanism  for  hydrogen  evolution  and 
oxidation.  The  mixed  potential  at  the  anode  is  calculated  using  the  mixed-potential  theory.  The  model 
results  are  compared  to  experimental  data  across  a  range  of  operating  conditions  and  component 
properties,  including  the  reactant  concentrations,  the  anolyte  flow  rate,  the  ionomer  volume  fractions 
and  the  membrane/ionomer  conductivity.  A  good  qualitative  fit  to  the  experimental  data  is  demon¬ 
strated.  In  order  to  gain  insight  into  the  anode  reaction  mechanism,  the  performance  on  both  Pt  and  Ni 
anodes  is  simulated  and  compared  to  experimental  observations.  A  detailed  analysis  of  the  results 
provides  an  explanation  for  the  different  performance  on  these  catalysts. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Direct  hydrogen  fuel  cells  have  the  potential  to  provide  clean 
and  sustainable  power  for  a  broad  range  of  applications,  from 
automotive  propulsion  to  portable  consumer  electronics  [1],  Cost 
effective  and  safe  methods  for  sourcing,  storing  and  handling 
hydrogen  are  key  to  the  mass  commericalisation  of  hydrogen  fuel 
cell  technology.  Widespread  adoption  also  requires  radical  changes 
to  existing  fuel  infrastructures.  In  order  to  alleviate  many  of  the 
problems  related  to  the  introduction  of  a  hydrogen  infrastructure, 
reformers  can  be  employed  to  convert  a  hydrocarbon  liquid  fuel  to 
hydrogen,  but  this  introduces  additional  engineering  complications 
and  increases  the  costs,  space  requirements  and  weights  of  the 
fuel-cell  power  systems. 
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Replacing  hydrogen  with  methanol  is  a  potential  solution  to  the 
aforementioned  problems.  Direct-methanol  fuel  cells  (DMFC),  on 
the  other  hand,  suffer  from  methanol  crossover  and  low  open- 
circuit  voltages  (OCVs)  [2],  in  part  due  to  the  inefficient  anodic 
oxidation  of  methanol.  The  problems  arising  from  the  use  of  liquid 
methanol  have  led  to  a  search  for  alternative  hydrogen-carrier 
compounds,  of  which  borohydride  compounds  are  a  prominent 
example  [3],  A  cell  that  utilizes  a  borohydride  compound  (usually 
sodium  borohydride)  as  the  fuel,  a  cation  or  anion  exchange 
membrane  and  02  [4]  or  H2O2  [5]  (in  anaerobic  environments)  as  an 
oxidant  is  termed  a  direct  borohydride  fuel  cell  (DBFC).  A  Schematic 
of  the  O2  DBFC  is  shown  in  Figure  1. 

DBFCs  are  characterized  by  their  high  OCVs  (considerably  higher 
OCVs  than  all  other  fuel  cell  technologies),  low  operating  temper¬ 
atures  and  high  power  densities  [6,7],  These  characteristics,  in 
combination  with  their  environmentally  friendly  by-products  of 
reaction,  have  generated  considerable  interest  in  DBFCs.  The  cell 
power  density  and  charge  efficiency  of  DBFCs  are  strongly 
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Fig.  1.  Schematic  of  a  direct  borohydride  fuel  cell  operating  with  02  as  the  oxidant. 


influenced  by  the  choice  of  anode  catalyst.  A  complete  (or  close  to 
complete)  eight  e  anodic  oxidation  is  feasible  with  Au  and  Ag 
electrodes  [3,8,9],  while  Pt  appears  to  liberate  up  to  only  four 
electrons  [10],  One  of  the  most  appealing  aspects  of  DBFCs  is  the 
ease  with  which  the  BH4  ions  can  be  oxidised  on  non-precious 
materials  such  as  Ni  and  Pd,  as  well  as  many  common  alloys  used 
for  hydrogen  storage  [7,10], 

Amongst  other  factors,  the  electrode  and  support  materials 
[7,10],  the  reactant  concentrations  [11,12],  the  PTFE/Nafion  binder 
content  in  the  anode  [13,14],  the  operating  conditions  [11,15],  the 
fuel  pH  [13],  the  diffusion  layer  properties  [16],  the  membrane  [7], 
the  electrode  structure  [17]  and  the  cell/stack  design  [7,10,18] 
influence  the  power  density  and  efficiency  of  DBFCs.  The  anode 
materials,  anolyte  composition,  membrane  and  operating  condi¬ 
tions  have  a  major  effect  on  the  evolution  of  H2,  which  is  an 
inherent  feature  of  DBFCs.  H2  evolution,  which  can  occur  via  at  least 
two  pathways,  lowers  the  charge  efficiency  of  the  cell,  inhibits 
reactant  transport  and  poses  safety  concerns.  The  cathode  perfor¬ 
mance  is  also  not  trivial.  It  has  been  demonstrated  in  several 
studies  that  the  cathode  polarisation  is  a  major  source  of  voltage 
loss  [4,12,19],  partly  due  to  the  establishment  of  a  mixed  potential 
resulting  from  borohydride  crossover  [20,21  ]. 

The  performance  of  current  DBFCs  is  highly  encouraging,  but 
further  improvement  and  optimization  are  required.  Amongst 
other  considerations,  there  is  ongoing  research  into  the  anode 
design,  particularly  with  respect  to  suppressing  H2  evolution, 
examining  and  understanding  the  voltage  losses  under  different 
design  and  operating  conditions,  and  further  understanding  the 
anode  and  cathode  reactions.  In  contrast  to  other  types  of  fuel  cells 
[22],  there  have  been  very  few  attempts  to  model  DBFCs  mathe¬ 
matically  [23],  Verma  and  Basu  [24]  developed  a  highly  simplified 
steady-state  model  of  an  02  DBFC  to  predict  the  cell  voltage  at 
a  given  current  density.  Activation  overpotentials  were  estimated 
by  inverting  Tafel  relationships  and  oxygen  diffusion  across  a  gas 
diffusion  layer  was  incorporated  by  linearising  Fick’s  law.  Diffusion 
in  the  catalyst  layer,  which  is  the  overwhelming  source  of  mass 
transfer  limitations,  was  neglected.  Ohmic  resistances  were 
modelled  by  estimating  an  overall  resistance  from  experimental 
data.  To  obtain  a  reasonable  qualitative  fit  to  experimental  data,  an 
ad  hoc  term  representing  the  concentration  polarisation  in  the 
anode  was  required.  Sanli  et  al.  [25]  developed  a  series  of  similar 
models  for  a  BH4  /H202  DBFC,  ignoring  mass,  charge  and  heat 
transport.  Concentration  overpotentials  were  derived  by  including 


reactant  concentrations  inside  a  Tafel  expression  and  by  intro¬ 
ducing  a  limiting  current  density.  The  ohmic  losses  were  charac¬ 
terised  by  a  constant  resistance,  estimated  experimentally  along 
with  the  reaction  constants  and  the  OCV.  Only  by  including  the 
cathode  explicitly  was  a  reasonable  fit  to  the  data  possible.  As  in 
Ref.  [24],  a  narrow  range  of  current  density  was  considered  (up  to 
only  0.02  A  cm-2).  Models  of  this  level  of  simplicity  can  provide 
a  means  for  monitoring  performance  but  lack  sufficient  physical 
detail  for  enhancing  fundamental  understanding  of  DBFC  operation 
and  for  improving  designs.  Material  properties,  solution  composi¬ 
tions,  the  complex  reaction  kinetics  and  the  operating  conditions 
are  known  to  exert  a  strong  influence  on  the  performance  of  DBFCs. 
Incorporating  these  aspects  of  DBFC  design  and  operation  is  key  to 
developing  effective  modelling  and  simulation  tools. 

In  this  paper,  a  physics-based  model  for  an  02  DBFC  is  devel¬ 
oped.  The  model  incorporates  the  main  factors  affecting  the  cell 
voltage.  Borohydride  oxidation,  hydrolysis  and  hydrogen  evolution/ 
oxidation  are  included  in  the  anode  mechanism.  The  full 
Tafel-Volmer-Heyrovsky  mechanism  is  employed  for  the 
hydrogen  evolution/oxidation  reactions,  together  with  competitive 
adsorption  of  H  and  BH4 .  The  mixed  potential  theory  is  employed 
to  estimate  the  anode  polarization.  The  predicted  trends  with 
respect  to  variations  in  key  variables,  including  the  reactant 
concentrations,  the  anolyte  flow  rate,  the  ionomer  volume  fractions 
and  the  membrane/ionomer  properties,  are  then  compared  to 
experimental  observations  (from  the  open  literature)  and  the 
results  are  analyzed  and  explained.  The  performance  on  a  Pt/C 
anode  and  the  performance  on  Ni  anode  are  also  simulated  and 
compared.  Explanations  for  the  differences  observed  in  the  simu¬ 
lations  (and  in  experiments)  are  proposed. 

2.  Mathematical  model 

Given  a  constant  applied  current  density  japplied.  the  output  of 
the  model  will  be  a  value  for  the  cell  voltage.  By  varying  jappiied. 
a  polarization  curve  can  be  simulated.  The  cell  voltage  is  calculated 
by  subtracting  the  activation,  ohmic  and  concentration  polar¬ 
isations  from  the  OCV  (at  which  the  net  current  flow  from  the 
anode  to  the  cathode  is  zero).  The  concentration  polarisations  are 
included  in  the  activation  polarisations  by  explicitly  including  the 
mass  transfer  resistances. 

In  the  model  developed,  the  cathode  is  considered  to  consist  of 
a  Pt/C  and  PTFE  catalyst  layer  pasted  onto  a  gas  diffusion  layer 
(GDL),  in  the  form  of  a  graphite  paper.  The  anode  consists  of  either 
Pt/C  or  Ni  mixed  with  PTFE  pasted  onto  a  Ni  foam.  The  electrodes 
are  separated  by  a  Nation  membrane.  A  schematic  of  the  MEA  is 
presented  in  Fig.  2.  The  cell  is  similar  to  that  developed  by  Li  et  al.  in 
Ref.  [  14]  and  where  parameter  values  are  available  from  that  paper, 
they  will  be  used  in  this  study.  It  is  important  to  point  out  that  the 
focus  of  this  study  is  to  develop  a  physics-based  framework  for 
modelling  DBFC  and  to  validate  the  approach  by  capturing  quali¬ 
tative  trends,  rather  than  quantitatively  fitting  to  experimental 
data.  There  is  no  loss  of  generality  in  this  approach;  once  suitable 
data  is  available,  together  with  all  material  properties,  a  quantita¬ 
tive  fit  would  be  straightforward  to  achieve. 

Several  simplifying  assumptions  were  adopted  in  the  develop¬ 
ment  of  the  model.  The  main  assumptions  and  their  justifications 
are  listed  below. 

1.  The  reactant  concentrations  are  assumed  to  be  spatially 
homogeneous  in  each  component,  although  the  main  transport 
limitations  are  incorporated. 

2.  The  generated  H2  gas  forms  gas  bubbles,  which  travel  at 
a  different  velocity  from  the  bulk  anolyte  velocity;  departures 
are  induced  by  drag  forces,  equal  in  magnitude  to  the  buoyancy 
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Bipolar  plate 


Fig.  2.  Schematic  of  the  cell  modelled  (the  components  enclosed  in  the  dotted  line). 


force  exerted  on  the  bubbles  [26],  Incorporation  of  these  effects 
is  not  attempted  in  this  paper.  The  H2  is  treated  as  either 
a  dissolved  or  an  adsorbed  species. 

3.  The  reaction  mechanism  in  the  anode  is  dependent  on  the 
catalysts  and  at  present  is  poorly  understood.  Therefore, 
a  general  mechanism  that  accounts  for  direction  borohydride 
oxidation,  hydrogen  evolution  and  hydrogen  oxidation  is 
adopted  (details  are  provided  in  the  next  subsection). 


2.1.  Reaction  mechanisms  and  reaction  rates 

Activation  losses  or  polarizations  are  a  consequence  of  the 
energy  barriers  associated  with  the  charge  transfer  reactions 
occurring  at  an  electrode.  In  other  words,  a  portion  of  the  voltage 
generated  is  consumed  in  transmitting  electrons  to  or  from  the 
electrode.  At  the  anode,  borohydride  ions  BHd  are  (ideally) 
oxidized  directly  on  the  electrode.  This  process  can  theoretically 
liberate  eight  electrons  [6,7]: 

BH4  +  80H  — ►  B02  +  6H20  +  8e  (e°  =  -1.24V  vs.  SHE)  (1) 

The  borohydride  fuel,  however,  tends  to  hydrolyze  rapidly  on  most 
candidate  electrocatalysts,  leading  to  the  formation  of  H2  [6,7]: 

BH4  +  2H20— >B02  +  4H2  (2) 

This  competitive  reaction  reduces  the  charge  efficiency  of  the  cell 
and  the  accompanying  formation  of  H2  leads  to  a  deterioration  in 
cell  performance  due  to  the  presence  of  H2  bubbles  in  the  anolyte 
(H2  molecules  hinder  ion  movement).  Reaction  (2)  is  accompanied 
by  a  second  pathway  to  H2  evolution,  namely  the  reduction  of 


Mechanism  (3)  is  the  full  Tafel— Volmer— Heyrovsky  mechanism,  in 
which  Hads  represents  an  adsorbed  adatom  of  hydrogen.  The  k ,  and 
k  t  are  forward  and  backward  reaction  rate  constants,  respectively. 
Water  reduction  is  thermodynamically  favoured  at  the  formal 
potential  for  borohydride  oxidation  (1),  suggesting  that  a  mixed 
potential  is  established  at  the  anode.  Indeed,  the  observed  anode 
potentials  are  significantly  more  positive  than  -1.24  V  vs.  SHE  and 
more  negative  than  -0.828  V  vs.  SHE  [7], 

There  is  evidence  that  Au,  Ag  and  some  hydrogen  storage  alloys 
are  capable  of  supporting  a  full  anodic  oxidation  of  the  borohydride 
under  certain  conditions  (either  due  to  inactivity  towards  hydro¬ 
lysis  or  possible  H2  adsorption  followed  by  oxidation),  while  Pt  and 
Ni  can  generally  support  only  a  four-electron  oxidation  [8,27—29], 
The  oxidation  of  certain  intermediates  of  hydrolysis,  such  as 
BH3(OH)“,  is  also  possible.  From  the  early  work  by  Bard  and  co¬ 
workers  [30,31],  the  initial  step  on  Au  is  thought  to  be  an  elec¬ 
tron  transfer  followed  by  a  rapid  decomposition  of  a  radical  and 
a  further  electron  transfer  step,  i.e.,  an 
electrochemical— chemical— electrochemical  sequence.  The  inter¬ 
mediates  (e.g.,  monoborane)  undergo  further  reactions  to  liberate 
a  six  (eight  in  total)  electrons,  though  the  precise  nature  of  these 
reactions  is  still  subject  to  speculation.  Chatenet  et  al.  [8]  studied 
the  direct  oxidation  of  NaBH4  on  Ag  and  Au  in  concentrated  NaOH 
solutions  using  linear  and  cyclic  voltammetry.  They  found  that  BH4 
electro-oxidation  depends  strongly  on  the  ratio  of  [BHd]  to  [OH-]. 
For  a  low  BH4  concentrations,  the  rate  of  chemical  hydrolysis  of 
BH4  was  negligible  but  at  low  ratios  of  [BH4  ]  to  [OH  ],  borohydride 
oxidation  was  limited  and  the  borohydride  ions  spontaneously 
hydrolyzed  to  BH3(OH)-,  which  undergoes  oxidation  at  low 
potentials.  There  are,  however,  a  number  of  apparently  contradic¬ 
tory  results  and  conclusions  in  the  literature,  in  part  due  to  the 
different  conditions  and  experimental  procedures  used  (such  as  full 
cell  vs.  rotating  disk  electrode).  For  example,  the  results  of  Li  et  al. 
[19]  suggest  that  full  anodic  oxidation  is  possible  on  Pt  under 
certain  conditions.  For  a  comprehensive  discussion  on  the  anode 
reactions  on  different  catalysts,  the  reader  is  referred  to  the  recent 
review  by  Merino-Jimenez  et  al.  [23]. 

Based  on  H2  evolution  data  for  a  range  of  anode  catalysts,  Wang 
et  al.  [32,33]  proposed  a  reaction  mechanism  consisting  of:  (a) 
a  four-electron  direct  borohydride  oxidation  involving  H2  evolution 
(reaction  (4)),  and  (b)  separate  hydrogen  evolution/oxidation  via 
reaction  (3): 

BH4,ads  +  40 H  ->•  B02  +  2H20  +  4e-  +  4Hads 

(h°  =  -  1 .24  V  vs.  SHE)  l  4) 

Reaction  (4)  is  a  special  case  of  the  reaction  proposed  by  Li  et  al.  [19] 
in  which  four  (of  a  possible  eight)  electrons  are  transferred  and  the 
H2  is  assumed  to  be  adsorbed.  While  this  mechanism  may  ignore 
the  possible  role  of  intermediates  such  as  BH3(OH)-,  it  provides 
a  framework  for  explaining  some  of  the  observed  effects,  such  as 
a  mixed  anode  potential.  The  model  developed  in  this  paper  is. 


fcv 

Volmer :  H20  +  e-  „  -  OH+Hads  i 

Vv  I 

Heyrovsky  :  Hads  +  H20  +  e-  J^V  H2  +  OH-  1  (e°  =  -0.828  V  vs.  SHE)  (3) 

Tafel :  2Hads  ==  H2  J 

kT 
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therefore,  based  on  the  mechanism  proposed  by  Wang  et  al. 
[32,33],  with  a  preliminary  step  of  BH4  adsorption,  as  proposed  by 
Chatenet  et  al.  [8]  (but  not  by  Wang  et  al.): 

BH4  ^  BH43ds  (5) 

VB„ 

The  rate  of  the  Tafel  step  in  mechanism  (3)  and  the  rate  of  reaction 
(5)  can  be  expressed  as  follows: 

rT  =  MptaH2  -  Mh^ 

rBH  =  ^  -  ^BH^BH^ 

respectively,  in  which  aBH-  and  aH2  are  the  surface  activities  of  BH4 
and  H2,  and  8m-,  and  dpt  are  the  catalyst-surface  coverages  of 
BH4 ,  the  H  adatoms  and  free  catalyst  sites,  respectively.  The  k  f  and 
k  j  are  again  the  forward  and  backward  rate  constants  for  the  two 
reactions.  As  a  simplification,  the  activity  (pressure)  of  hydrogen  in 
equation  (6)  is  assumed  to  be  unity. 

The  Butler— Volmer  equation  [34]  can  be  used  to  approximate 
the  transfer  current  density  associated  with  reaction  (4).  Cheng  and 
Scott  showed  that  the  transfer  current  density  is  first-order  with 
respect  to  OH-  [35],  Since  the  concentration  of  water  in  the  anolyte 
is  much  higher  than  the  concentration  of  borohydride  ions  and 
sodium  hydroxide,  it  is  reasonable  to  assume  a  unit  activity  for 
water.  The  dependence  on  the  B02  concentration  is  also  neglected 
since  the  reverse  of  reaction  (4)  is  exponentially  small  at  the  mixed 
anode  potential.  Accordingly: 

-exp(-('-y»)}  (7) 

where  j'bh  is  the  borohydride  oxidation  current  density;  Jbh, ref  is  the 
reference  exchange  current  density  associated  with  the  borohy¬ 
dride  oxidation  reaction;  c0H-  is  the  concentration  of  OH-  at  the 
electrode  catalyst  surface  and  cgj£-  is  the  corresponding  reference 
concentration  (at  which  j'bh  —  jBH.ref);  Vm  is  the  overpotential;  R  is 
the  universal  gas  constant;  T  is  the  cell  temperature;  Fis  Faraday’s 
constant;  na  is  the  number  of  electrons  transferred  (four);  and  a\  is 
the  charge  transfer  coefficient.  The  overpotential  is  given  by: 

^BH  =  0s  “  0e  “  EfiH  (8) 

in  which  is  the  equilibrium  potential  for  reaction  (4),  <j>s  is  the 
electronic  potential  and  <j>e  is  the  ionic  potential. 

The  transfer  current  density  associated  with  the  water  reduction 
reaction  (3)  is  estimated  using  a  generalized  form  of  the 
Butler— Volmer  equations  for  the  Volmer  and  Heyrovsky  steps 
[36—40],  assuming  a  unit  water  activity  and  assuming  that  the 
sodium  hydroxide  is  not  limiting: 

Jv  =  F{^y°»exp(a2RTH)  ~  fevflptexp^1  ~g°2)F%) ] 

(9) 


(10) 


where  77H  is  the  overpotential  and  a2  is  the  charge  transfer  coeffi¬ 
cient.  A  unit  hydrogen  activity  is  again  assumed  in  equation  (10). 
The  overpotential  is  defined  as: 

?1H  =  0S-0e-£Hq  0D 

where  E^q  is  the  equilibrium  potential  for  reaction  (3). 

The  reaction  in  the  air-breathing  cathode  is  represented  as 
follows  [6,7]: 

02  +  2H20  +  4e~  — >40H~  (e°  =  0.4  V  vs.  SHE)  (12) 

The  cathode  transfer  current  density  can  be  expressed  as 
follows,  making  the  reasonable  approximation  of  a  unit  water 
activity: 

fa  =w(|){exp(!^)-exp(=M^)} 

(13) 


where  jo2  ref  is  the  exchange  current  density  associated  with  reac¬ 
tion  (12);  Co2  is  the  oxygen  concentration  at  the  reaction  sites 
(platinum  surfaces)  and  cge2f  is  the  corresponding  reference 
concentration;  nc  is  the  number  of  electrons  transferred;  and  rjc  is 
the  cathode  overpotential.  The  dependence  on  OH-  is  neglected 
since  reaction  (12)  is  irreversible  and  the  OH-  ions  are  not  limiting. 
The  cathode  overpotential  is  defined  as: 

Vc  =  0s-0e-£cq  (14) 

where  E^q  is  the  equilibrium  potential  for  reaction  (12). 

2.2.  Activation  overpotentials  and  anode  polarization 


For  a  non-zero  constant  applied  current  density,  jappiied,  the 
cathode  overpotential  is  obtained  by  inverting  the  following 
equation,  which  relates  the  transfer  current  density,  j02,  defined  in 
equation  (13),  to  jappiied: 

rd02  =  Jappiied  (15) 


where  the  quantity  rc  =  A“th/Ag  is  the  roughness  factor  of  the 
cathode;  A“th  is  the  active  surface  area  for  reaction  in  the  cathode 
and  Ag  is  the  geometrical  surface  area  of  both  electrodes  (the 
exchange  current  density  is  per  unit  area  of  the  active  catalyst, 
while  the  applied  current  density  is  per  unit  geometric  area  of  the 
electrode).  An  expression  for  the  surface  oxygen  concentration  is 
derived  in  the  next  subsection. 

The  potential  at  the  anode  is  a  mixed  potential  [41  ]  due  to  the 
simultaneous  occurrence  of  the  borohydride  oxidation  and  water 
reduction  reactions  (4)  and  (3).  Under  the  condition  of  a  zero 
current  density,  these  reactions  are  not  at  equilibrium  since  the 
mixed  potential  E°ixed  lies  between  the  equilibrium  potentials  for 
reactions  (4)  and  (3),  with  the  anodic  branch  of  (4)  and  the  cathodic 
branch  of  (3)  dominating,  i.e.,  there  is  a  balance  between  the 
electrons  released  in  reaction  (4)  and  electrons  consumed  in 
reaction  (3).  Mathematically,  this  can  be  expressed  as  follows,  using 
equations  (7),  (9)  and  (10)  with  a^  =  a2  =  0.5  [41]: 


Dr 


n3F(E°mx 


-F("kv0pt  +  DH«H)exp 


^(£mixed  —  iff) 


Jh  =  F | Dh 0ptexp -7HtfHexp(  !  l  r“^F7?H)} 


(16) 
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Given  values  for  the  reactant  concentrations  and  for  £j^  and  £^q, 
equation  (16)  can  be  solved  for  the  mixed  anode  open-circuit 
potential  EjF^  (i.e.,  the  difference  between  the  electronic  and 
ionic  potentials,  </>s  -  <t>e,  at  zero  current  density).  For  a  non-zero 
constant  applied  current  density,  jappiied.  conservation  of  elec¬ 
tronic  charge  demands  that: 

ra(/BH  +JV  +Jh)  =  Jappiied  C17) 

where  the  quantity  ra  =  Af'/Ag  is  the  roughness  factor  of  the 
anode.  Therefore: 

+  k  v$HexP  (^2pf  (^mixed  — 

-  Tvf?Ptexp(  -  2Frt  (£mixed  -  £*q))  ]  (18) 

+  f{  tnAptexp  (£mixed  -  O  ) 

—  k  H^HexP^  —  2pj  (^mixed  —  J" 

Equation  (18)  can  be  solved  for  the  mixed  anode  potential  £mixed. 
which  is  the  difference  between  the  electronic  and  ionic  potentials 
pertaining  to  a  non-zero  current  density.  The  polarization  of  the 
anode,  i]mjxed  is  calculated  from  the  following  equation: 

Wd  =  £mixed- C<ed  (19) 


2.3.  Mass  balances 


Calculation  of  the  activation  losses  in  both  the  anode  and 
cathode  requires  estimates  of  the  reactant  concentrations  at  the 
electrode  surfaces,  where  the  reactions  take  place.  The  DBFC  in  Ref. 
[4]  utilized  porous  electrodes  for  both  the  anode  and  cathode.  In 
the  model  developed,  the  anolyte  is  circulated  through  the  anode  at 
a  constant  volumetric  flow  rate  w.  A  liquid  diffusion  layer  (Ni  foam) 
is  assumed  to  separate  the  anode  catalyst  layer  (ACL)  and  the 
anolyte  inlet  and  outlet  (a  common  design).  It  is  assumed  that  the 
electric  current  enters  uniformly  along  bipolar  plate  surfaces,  i.e.,  in 
the  through-plane  direction  x  in  Fig.  2. 

The  concentrations  of  the  active  species  in  the  anode 
(comprising  the  catalyst  layer  and  Ni  foam)  are  derived  from  a  mass 
balance  based  on  the  circulation  of  the  anolyte  together  with 
reactant  consumption  or  generation.  At  practical  flow  rates,  diffu¬ 
sive  transport  in  the  anode  is  negligible  in  comparison  to  convec¬ 
tive  transport,  except  within  mass-transfer  boundary  layers.  To 
model  the  resistance  to  mass-transfer  between  the  bulk  region 
(pore  space)  and  the  reaction  sites  (catalyst  surfaces),  a  Nernst 
diffusion  boundary  layer  [41]  is  introduced.  The  bulk  anolyte 
solution  is  considered  well  mixed  in  the  catalyst  layer  and  Ni  foam, 
leading  to  uniform  bulk  species  concentration  values  that  are  valid 
in  both  regions.  At  steady  state,  molar  balances  for  the  species  i  = 
BHJ ,  BO2  ,  OH-  in  the  bulk  region  can  be  written  as  follows: 


fc“k  -  c, 


(20) 


/  Met  bulk!  AlnDf  ( 

w(c!  -CH  -  V 

where  A“  is  the  surface  area  to  which  the  reactants  diffuse  through 
the  boundary  layer  (the  active  surface  area  for  reaction  in  the 
anode);  <5,-  is  the  Nemst  boundary  layer  thickness  for  species  i;  c[nlet, 
cbu!k  ancj  c.  are  the  inlet,  bulk  (in  the  catalyst  layer  and  Ni  foam)  and 
catalyst  surface  concentrations  of  species  i,  respectively;  and  D?ff  is 


the  effective  diffusion  coefficient  of  species  t  in  solution  in  the 
catalyst  layer.  Dfff  can  be  approximated  using  the  Bruggeman 
correction  [42]:  Dfff  =  where  raci  is  the  porosity  of  the  ACL 

and  Dj  is  the  free-space  diffusion  coefficient  in  solution.  The  steady 
state  molar  balances  for  species  i  —  BO2  ,  OH-  in  the  diffuse  layer 
can  be  written  as  follows: 

vpf  (c-  _  Cl)  _  .0  (2,) 


in  which  r,-  is  the  stoichiometric  coefficient  of  species  i  in  reaction 
(4),  and  pi  and  tt,-  are  the  stoichiometric  coefficients  of  species  i  in 
the  Volmer  and  Heyrovsky  steps  (3),  respectively  (taken  as  positive 
if  the  species  is  a  reactant  and  negative  if  the  species  is  a  product). 
Eliminating  the  bulk  concentration  in  equations  (20)  and  (21) 
yields  expressions  for  the  surface  and  bulk  concentrations  of 
species  i  —  BO2 ,  OH- : 


inlet  ( AT  ,  Mm  Wv  +  fftflA 

'■  +DfffJW  F  ) 


The  surface  concentration  of  BH4  is  calculated  from  a  steady- 
state  balance  that  takes  into  account  diffusion  to  the  reaction 
sites  and  adsorption  and  desorption,  the  rates  of  which  are  given  in 
equation  (6).  The  activity  of  BH4  in  the  latter  equation  is  approxi¬ 
mated  by  the  ratio  of  the  BH4  surface  concentration  to  the  inlet 
concentration,  c^et.  The  mass  balance  therefore  reads: 


-  (cBH„k  “  cBH4)  =  ^anfBH 

=  Aan(kBH^PtCBH4  -  /<BH0BH4) 


which  is  combined  with  equation  (20).  In  equation  (23),  the 
modified  reaction  constant  kBH  =  k  eh/Cebt  has  been  introduced 
to  simplify  notation.  Equations  are  also  required  for  the  surface 
coverages  of  hydrogen  and  BH4 ,  which,  under  steady  state  condi¬ 
tions,  are  given  by: 


(24) 


respectively,  together  with: 

0pt  +  0H  +  0BH4  =  1  (25) 


i.e.,  conservation  of  the  active  catalyst  sites. 

There  are  various  approaches  to  modelling  the  diffusion 
boundary  layer  thickness  d.  For  each  oxidized/reduced  species  in  an 
electrolyte— electrode  structure,  Chrysikopoulos  et  al.  developed  an 
expression  for  the  boundary  layer  thickness  in  terms  of  the  (mean) 
flow  rate  in  a  direction  parallel  to  the  layer  surface  [43]: 


sfi/Df 


(26) 


where  K  is  a  constant  parameter  depending  on  the  area.  The 
formula  (26)  is  adopted  in  this  paper. 

The  sodium  ions  do  not  participate  in  the  anode  reactions.  They 
do,  however,  migrate  from  the  anode  to  the  cathode  under  the 
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influence  of  the  potential  gradient  formed  under  a  load.  Assuming 
that  BHJ,  BO2  ■  OH-  and  Na+  are  the  only  charged  species  in  the 
electrolyte,  the  concentration  of  Na+  can  be  calculated  from  the 
condition  of  electroneutrality  in  the  bulk  solution  (close  to  the 
electrode  surface,  this  condition  is  of  course  violated  and  a  double 
layer  develops): 

J2z4u,k  =  0  (27) 


where  zi  is  the  charge  number  of  species  i.  Equation  (27)  ensures 
that  (in  the  model)  Na+  migrates  to  the  cathode  through  the 
membrane  in  order  to  maintain  charge  balance  in  the 
electrodes. 

The  oxygen  is  fed  to  the  cell  through  a  cathode  flow  channel.  It 
diffuses  through  a  gas  diffusion  layer  (GDL)  to  the  cathode  catalyst 
layer  (CCL)  where  it  reacts  on  the  Pt/C  catalyst.  In  such  air-breathing 
Pt/C  designs,  the  O2  must  diffuse  through  the  ionomer  in  the 
catalyst  layer  in  order  to  reach  the  reaction  sites,  which  introduces 
an  additional  mass-transport  resistance.  This  resistance  can  be 
incorporated  by  defining  a  boundary  layer  of  thickness  50l .  A  steady 
state  molar  balance  of  oxygen  in  the  catalyst  can  be  derived  by 
taking  into  account  O2  diffusion  through  the  GDL  to  the  CCL,  where 
the  concentration  is  denoted  c^ulk,  mass  transfer  through  the 
boundary  layer  to  the  catalyst  surfaces,  where  the  concentration  is 
denoted  c02,  and  electrochemical  reaction: 


pcff/j  C°“' 

°o/g  /cdT~ 

^cathpeff  k  \  Japplied^g  „ 

<5n(  v  O2  -Co2J - 4F  =  0 


AcathDeff 

"V-(C“-O  =  o 


In  equation  (28):  A“th  is  the  active  surface  area  for  reaction  in  the 
cathode  (the  active  catalyst  surface  area,  to  which  the  reactants 
diffuse  through  the  boundary  layer);  Igdl  is  the  thickness  of  the 
GDL;  Dfff  is  the  effective  diffusion  coefficient  of  oxygen  through  the 
ionomer  film  (boundary  layer)  and  is  the  effective  diffusion 
coefficient  of  oxygen  through  the  GDL.  The  effective  diffusion 
coefficients  are  approximated  using  Bruggeman  corrections  [42]: 
Df  =  Of  and  Df2  =  42lD02,  in  which  £e,cci  is  the  volume 
fraction  of  ionomer  in  the  CCL,  £Gdl  is  the  porosity  of  the  GDL,  Df  is 
the  diffusion  coefficient  of  oxygen  in  the  bulk  ionomer  and  Dq2  is 
the  free  space  diffusion  coefficient  of  oxygen  in  air.  Eliminating  the 
bulk  concentration  in  equation  (28)  yields  an  expression  for  the 
surface  concentration  of  oxygen: 


co2 


*GDL  .  ^02  \  ./applied 

rcDf  J  4 F 


(29) 


2.4.  Ohmic  losses 

The  resistances  to  the  flow  of  charge  in  the  ion-conducting  and 
electron-conducting  phases  result  in  (additive)  voltage  losses, 
which  are  calculated  using  Ohm’s  law.  The  relevant  regions  are  the 
ACL,  the  Ni  foam,  the  membrane,  the  CCL  and  the  cathode  GDL.  Both 
the  liquid  and  solid  phases  are  taken  into  account. 

2.4.1.  Ionic  resistances 

The  ohmic  drop  across  the  membrane  can  be  written  in 
terms  of  the  Na+  conductivity  of  the  membrane,  <7mem,  and  the 
applied  current  density  (the  current  is  purely  ionic  in  the 
membrane): 


(30) 


in  which  /mem  is  the  thickness  of  the  membrane. 

Ion  transport  through  the  ionomer  in  the  ACL  (it  is  assumed  that 
the  ionomer  is  the  same  material,  e.g.  Nation,  as  the  membrane 
material,  which  is  typically  the  case)  takes  place  alongside  ion 
transport  through  the  anolyte  solution.  The  ionic  potential  drop 
across  the  ACL  must  take  into  account  these  parallel  processes.  The 
drop  in  the  ionic  potential  across  the  ACL  can  be  calculated  by 
considering  the  total  ionic  current  density  in  the  in  the  ACL,  je, 
generated  by  the  flow  of  charged  species  in  the  anolyte  and  the 
migration  of  sodium  ions  in  the  ionomer  (in  the  through  plane 
direction  x  in  Fig.  2)  [34]: 


In  this  equation,  v  is  the  electrolyte  flow  speed  and  N;  is  the 
molar  flux  of  species  i,  which  is  expressed  using  the  Nemst— Planck 
equation.  The  effective  diffusion  coefficients  in  the  ACL,  D?ff,  are 
approximated  using  the  Bruggeman  correction  [42]:  Dfff  =  e^D;. 
The  conductivity  of  the  ionomer  is  also  based  on  a  Bruggeman 
correction,  using  the  volume  fraction  of  ionomer  in  the  ACL,  Ee,acl- 
The  term  5DZjFt/cj,ulk  in  equation  (32)  is  zero  by  virtue  of  the  elec¬ 
troneutrality  condition  (27)  and  the  term  X)ZiFD?ff(dckulk/dx)  is 
zero  given  the  assumption  of  uniform  bulk  concentrations.  The 
term  rf^F2  ^z?DiC^u,k/(RT)  is  the  effective  conductivity  of  the 
anolyte  and  the  overall  effective  ionic  conductivity  in  the  ACL,  ff^cl 
is  given  by  the  sum  of  the  effective  conductivity  of  the  anolyte  and 
the  effective  conductivity  of  the  ionomer,  reflecting  the  parallel 
paths  for  ion  conduction.  Thus,  the  total  ionic  potential  drop  across 
the  ACL  is  given  by: 

(^e)ad  =  Jg-je  (32) 

°e,ad 

where  lac  1  is  the  thickness  of  the  ACL. 

For  the  ohmic  drop  across  the  CCL  it  is  again  assumed  that  the 
ionomer  is  the  same  material  as  the  membrane  material.  The  ohmic 
drop  can  then  be  written  in  terms  of  the  Na+  conductivity  of  the 
membrane  and  the  ionic  current  density: 

(A0e)ccl  =  3/2^  Je  (33) 


in  which  (cci  is  the  thickness  of  the  CCL.  Note  that  the  effective 
conductivity  has  been  approximated  by  ^  cciamem>  using  a  Brugge¬ 
man  correction  [42], 

2.4.2.  Electronic  resistances 

To  find  the  electronic  potential  drop  in  the  ACL,  the  volume 
fractions  of  catalyst  (either  Pt  or  Ni  in  the  simulations)  and  carbon, 
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t'M.a  and  tc.aci.  respectively,  and  the  individual  conductivities  of  the 
catalyst  and  carbon,  <tm  and  at,  respectively,  are  used: 


”  s  arl  r,i,aM  +  Cr- 


(34) 


where  the  effective  conductivity  a^cl  is  given  by  the  sum  of  the 
individual  conductivities  (with  Bruggeman  corrections),  corre¬ 
sponding  to  resistances  in  parallel.  The  current  density  through  the 
solid  electron-conducting  components  is  denoted  js  in  equation 
(35). 

The  electronic  potential  drop  across  the  CCL  can  be  derived  in 
a  manner  similar  to  that  leading  to  equation  (35): 


(A0s)ccl  ~  ofJs  ~  25s 


(35) 


where  the  effective  conductivity  a^cl  is  again  given  by  the  sum  of 
the  individual  conductivities  (of  carbon  support,  with  volume 
fraction  tc.cd.  and  platinum,  with  volume  fraction  tptiC  and 
conductivity  apt). 

The  electronic  potential  drop  across  the  Ni  foam  layer  (in  which 
the  current  is  purely  electronic)  is  given  by: 


(A0e)Ni 


(1  -  ENi)3/2(7Ni 


(36) 


where  In:  is  the  thickness  of  the  Ni  foam  and  aw  is  the  conductivity 
of  nickel.  The  electronic  potential  drop  across  the  (carbon)  cathode 
GDL  is  given  by: 


(A0s)gdl  -  — 


^GDL _ 

£‘gdl)3/2(7c 


(37) 


The  current  is  again  purely  electronic  through  this  component. 
The  charge  balances  are  not  explicitly  incorporated.  The  values  of  js 
and  je  are  instead  approximated  by  their  (formal)  averages,  that  is, 
js  =js  =  jappiied/2  in  both  catalyst  layers.  Further  details  are  provided 
in  the  Appendix. 


2.5.  Calculation  of  the  cell  voltage 

The  cell  voltage  can  be  calculated  by  subtracting  the  ohmic 
losses  and  the  activation  and  concentration  polarizations  from  the 
OCV: 

^cell  —  ^cell  3“  ^mixed  ~  Vc~  (^0e)ad  —  (^</’e)ccl  —  Welmem 
-  (A*,)*  -  (Ms)ccl  -  (^s)Ni  -  (^s)GDL 

(38) 

in  which  all  terms  except  E°eH  have  been  defined  in  Sections 
2.1 -2.4.  Note  that  the  concentration  (mass  transport)  polarization 
is  incorporated  in  the  polarizations  ijmixed  and  r]c,  which  are 
dependent  on  the  concentrations  of  the  reactants  at  the  electrode 
surfaces,  which,  in  turn,  incorporate  the  major  mass  transport 
resistances. 

The  OCV  is  the  difference  between  the  cathode  equilibrium 
potential  and  the  mixed  anode  open  circuit  potential  (Section  2.2): 

Ec°e„  =  Ef  -  E°  ixed  (39) 

The  individual  equilibrium  potentials  for  the  reactions  (4),  (3)  and 
(12)  will  be  required  to  evaluate  Ej?ell;  they  can  be  expressed 


Table  1 

Default  electrochemical  parameters  used  in  the  model. 


Symbol 

Quantity 

Value 

Reference 

E? 

Standard  potential  of 
oxygen  reduction 

0.4  V  vs.  SHE 

[7] 

E°H 

Standard  potential  of 
borohydride  oxidation 

-1.24  V  vs.  SHE 

[7] 

Eh 

Standard  potential  of 
hydrogen  oxidation 

-0.828  V  vs.  SHE 

[7] 

at 

Charge  transfer  coefficient 

0.5 

Assumed4 

°2 

Charge  transfer  coefficients 
for  reactions  (3) 

0.5 

[40] 

a3 

Charge  transfer  coefficient 

0.33 

[46] 

Jm,  ref 

Reference  exchange  current 
density  for  reaction  (4) 

60  A  m"z 

[45] 

% 

Backward  rate  constant  for 
Volmer  reaction  (3) 

4.4  x  10  3  mol  m2 

s-1  [40] 

T(v 

Forward  rate  constant  for 
Volmer  reaction  (3) 

4.4  x  10-4  mol  m2 

s-1  [40] 

Backward  rate  constant 
for  Heyrovsky  reaction  (3) 

2.4  x  10-6  mol  m2 

s-1  [40] 

t  H 

Forward  rate  constant  for 
Fieyrovsky  reaction  (3) 

2.4  x  10-5  mol  m2 

s-1  [40] 

tT 

tT 

Backward  rate  constant 
for  Tafel  reaction  (3) 

8.8  x  10-6  mol  m2 

s-1  [40] 

Forward  rate  constant 
for  Tafel  reaction  (3) 

8.8  x  10  4  mol  m2 

s-1  [40] 

*kBH 

Backward  rate  constant 
for  reaction  (5) 

0  mol  m2  s-1 

[8] 

~kBn 

Forward  rate  constant 
for  reaction  (5) 

1  x  10“5  mol  m2  s' 

Fitted 

JOt.'ef 

Exchange  current  density 
for  reaction  (12) 

2.1  x  10  3  A  m  2 

[46] 

r* 

Roughness  factor  of 

100 

Assumed4 

rc 

Roughness  factor  of  the 
cathode 

77 

[57] 

mr« 

Pt  loading  in  the  cathode 
catalyst  layer 

4gm-2 

Assumed4 

mew 

catalyst  layer 

40  g  m  2 

Assumed4 

4  See  Section  2.6. 


using  the  Nernst  equation  [41  ]  but  as  a  first  approximation  they  are 
taken  to  be  equal  to  their  standard  values  (Table  1 ). 

2.6.  Model  parameters 

The  default  parameter  values  for  the  model  are  shown  in 
Tables  1—4.  The  values  of  the  properties  labelled  ‘Assumed’  are 
based  on  the  materials  described  below,  which  are  typical  for 
DBFCs. 

The  cathode  GDL  properties  are  based  on  Toray  TGP-FI-060 
paper  and  the  Ni  foam  properties  are  based  on  a  commercially 
available  foam  from  Sumitomo  Electric.  The  anode  catalyst 


Table  2 

Default  operating  conditions  used  in  the  model. 


cirdel 

c£f 


Quantity 

Operating  temperature 
Inlet  concentration  of  BHJ 
Inlet  concentration  of  BOJ 
Inlet  concentration  of  OH~ 
Inlet  concentration  of  Na+ 
Inlet  concentration  of  02 
Reference  02  concentration 


298  K 
0.5  M 
0.0  M 
2.0  M 
2.5  M 

8.53  mol  nr3 
8.53  mol  m-3 
10  mL  min1 
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Table  3 

Default  structural  parameters  used  in  the  model. 


Symbol  Quantity 


Reference 


^Oj 

A. 


Porosity  of  cathode  GDL  (Toray) 

Ionomer  volume  fraction  in  cathode 
catalyst  layer 

Ionomer  volume  fraction  in  anode 
catalyst  layer 

Ni  foam  porosity  (Sumitomo  Electric) 

Volume  fraction  of  Pt  in  cathode  catalyst  layer 
Volume  fraction  of  Pt  in  anode  catalyst  layer 
Volume  fraction  of  carbon  in  the  cathode 
catalyst  layer 

Volume  fraction  of  carbon  in  the  anode 
catalyst  layer 

Thickness  of  cathode  GDL  (Toray) 

Thickness  of  Ni  foam  (Sumitomo  Electric) 
Thickness  of  cathode  catalyst  layer 
Thickness  of  anode  catalyst  layer 
Thickness  of  membrane  (Nafion®  112/212) 
Diffusion  film  thickness  in  the  cathode 
Geometric  surface  area  of  the  electrodes 


0.73  [58] 

0.2  [59] 

0.2  Assumed3 

0.98  Assumed3 
0.09  Assumed3 
0.09  Assumed3 
0.25  Assumed3 

0.31  Assumed3 

190  pm  [58] 

1.5  mm  [58] 

20  pm  [57] 

200  pm  Assumed3 
50  pm  Assumed3 
2  pm  Assumed3 
10  cm2  Assumed3 


3  See  Section  2.6. 


properties  in  the  tables  are  based  on  Pt/C.  Values  based  on  Ni  are 
provided  in  Table  5.  The  film  diffusion  coefficient  in  the  cathode,  Df 
(Table  4),  is  based  on  the  diffusion  coefficient  through  a  Nafion 
membrane  at  room  temperature.  The  reference  exchange  current 
density  for  borohydride  anodic  oxidation  on  Pt/C  (/BH.ref  in  equation 
(7))  was  calculated  at  25  °C  by  Molina-Concha  and  Chatenet 
[44,45],  ^ 

The  reaction  rate  constants  k ,  and  k ;  for  the  hydrogen  evolu¬ 
tion  reaction  (3)  on  Pt  (Table  1)  are  taken  from  Bai  [38]  (see  also 
[39,40]),  estimated  at  23  °C.  Equivalent  values  for  Ni  (Table  5)  were 
taken  from  the  work  of  Jaksic  et  al.  [37],  The  exchange  current 
density  and  transfer  coefficient  for  oxygen  reduction  (Table  1 )  were 
taken  from  Yu  et  al.  [46],  who  measured  the  values  for  60  °C  and 
25  °C.  The  reference  concentration  of  oxygen  was  assumed  to  be 
equal  to  the  inlet  value,  which  was  based  on  air  at  1  bar  and  25  °C  in 
the  default  case  (Table  2). 

The  volume  fractions  of  catalyst  in  the  anode  and  CCLs  (Table  3) 
can  be  found  from: 


P  M^ad 


mPt,c 

/Wccl 


(40) 


respectively,  where  him, a  and  mi>t,c  are  the  anode  catalyst  and 
cathode  Pt  loadings  (kg  per  m2  of  the  anode  geometric  surface 
area),  respectively,  pu  and  pn  are  the  densities  of  the  anode  catalyst 


Table  4 

Default  mass-  and  charge-transport  properties  used  in  the  model. 


Symbol 

Quantity 

Value 

Reference 

£>bh„ 

Diffusion  coefficient  of  BHJ 
in  anolyte  (2  M  NaOH) 

1.0  x  to  5  cm2 

s *  1  [28] 

D  BOJ 

Diffusion  coefficient  of  BOJ 
in  anolyte 

2.24  x  10  5  cm' 

!  s_1  [29] 

t>OH- 

Diffusion  coefficient  of  OH~ 
in  anolyte  (2  M  NaOH) 

8.2  x  10  6  cm2 

s_1  [60] 

DNa+ 

Diffusion  coefficient  of  Na+ 
in  anolyte  (2  M  NaOH) 

5.5  x  to  6  cm2 

s-1  [60] 

Doa 

Diffusion  coefficient  of  O2 
in  GDL  (in  air  @  STP) 

0.21  cm2  s— 1 

[29] 

Df 

Diffusion  coefficient  of  O2 
through  diffusion  film 
layer  (Nafion)  in  the  cathode 

5.9  x  10  7  cm2 

s'1  [61] 

"GDI. 

Effective  conductivity  of 
cathode  GDL 

12.5  S  cm1 

[58] 

K 

Constant  in  equation  (26) 

0.5  mV2 

Assumed3 

3  See  Section  2.6. 


Symbol  Quantity 

fli  Charge  transfer  coefficient 

for  reaction  (4) 

a2  Charge  transfer  coefficients 

for  reactions  (3) 

jBH.ref  Reference  exchange  current 

density  for  reaction  (4) 

for  Volmer  reaction  (3) 
fcv  Forward  rate  constant 

for  Volmer  reaction  (3) 
k  H  Backward  rate  constant  for 

Heyrovsky  reaction  (3) 
k  h  Forward  rate  constant  for 

Heyrovsky  reaction  (3) 
fcT  Backward  rate  constant 

for  Tafel  reaction  (3) 
kT  Forward  rate  constant 

for  Tafel  reaction  (3) 
k  bh  Backward  rate  constant 

for  reaction  (5) 

k  BH  Forward  rate  constant 

for  reaction  (5) 

ra  Roughness  factor  of 

mNi  a  Ni  loading  in  the  anode 

catalyst  layer 


6  x  10  5  mol  m2  s  1 
5  x  10~6  mol  m2  s'1 
4.2  x  10  s  mol  m2  s_1 
1.5  x  10~7  mol  m2  s_1 
7.8  x  10~8  mol  m2  s-1 
13  x  10~5  mol  m2  s-' 


1  x  10~6  mol  m2  s"1 


40  g  m~2 


3  See  Section  2.6 


Reference 

Assumed3 


[40] 

[45] 

[37] 

[37] 

[37] 

[37] 

[37] 

[37] 


Assumed3 

Assumed3 


and  cathode  Pt,  respectively.  The  default  loadings  in  Table  1,  both 
based  on  Pt,  were  taken  to  be  0.4  mg  cm-2  and  4  mg  cm2  for  the 
cathode  and  anode,  respectively,  both  with  40%  catalyst-to-carbon 
wt.%  (these  values  are  typical  of  unit  DBFCs). 

The  adjustable  (fitted)  parameters  are  listed  below. 

1.  The  reaction  constants  k  BH  and  k  Bh  in  equation  (6).  Following 
Chatenet  et  al.  [8],  the  backward  constant,  representing 
desorption  is  considered  zero,  i.e.,  adsorption  is  irreversible. 
The  selected  value  of  k  BH  (Table  1 )  will  be  discussed  in  the  next 
section. 

2.  The  boundary  layer  thickness  in  the  anode  is  specified  by  the 
value  of  the  constant  K  in  equation  (26).  The  value  of  I<  is 
system  specific  and  is  not  available.  The  selected  value  of  K 
(Table  4)  will  again  be  discussed  in  the  next  section. 

3.  The  Na+  conductivity  of  the  Nafion®  type  membrane  depends 
on  the  membrane  water  content  and  the  temperature  but  the 
precise  nature  of  the  dependence  has  not  been  characterized. 

The  conductivity  also  depends  on  the  concentration  of  the 
NaOH  in  the  anolyte  to  which  the  membrane  is  exposed  on  the 
anode  side.  The  results  of  Nargbska  et  al.  [42]  suggest  a  Nafion® 
120  conductivity  value  of  between  1  and  2.5  S  m-1  for  NaOH 
concentrations  in  the  range  0.1—2  M  in  aqueous  solution  at 
298  K.  Ma  et  al.  [47]  measured  the  conductivity  of  a  Nafion® 
212  membrane  equilibrated  with  10  wt.%  aqueous  NaOH 
solution  at  room  temperature  and  recorded  a  value  of 
0.74  S  m1.  A  value  of  0.5  S  nr1  (Table  4)  was  chosen  as  the 
default  value  and  the  effects  of  deviations  from  this  value  were 
investigated,  as  detailed  in  the  next  section. 


3.  Results  and  discussion 

3.1.  Performance  using  a  platinum  anode  catalyst 

Simulations  were  performed  using  the  default  values  in 

Tables  1—4  and  the  results  are  depicted  in  Figs.  3—5  (the  term 
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Current  Density  /  A  cm  2 

Fig.  3.  The  cell  voltage  and  the  individual  component  of  the  total  voltage  loss  vs.  the 
applied  (total)  current  density  for  the  base-case  parameters  in  Tables  1-4.  The  anode 
catalyst  is  Pt/C. 


‘Current  density’  on  the  x-axis  of  all  figures  refers  to  the  applied 
(total)  current  density).  Fig.  3  shows  the  cell  voltage  at  different 
applied  current  densities,  along  with  the  individual  components  of 
the  voltage  loss.  This  figure  clearly  indicates  that  the  ohmic  losses 
represent  a  significant  portion  of  the  total  losses,  which  is  consistent 
with  experimental  observations  for  both  O2  and  H2O2  systems,  e.g., 
[48—50],  The  ohmic  losses  are  dominated  by  the  resistance  to  ion 
transport  in  the  membrane,  the  ionomers  and  the  electrolyte  solu¬ 
tion,  as  discussed  later.  Another  notable  feature  of  Fig.  3  is  the 
relatively  small  anode  polarization  (deviation  of  the  mixed  anode 
potential  from  the  anode  mixed  open-circuit  potential),  particularly 
in  comparison  to  the  cathode  activation/mass-transport  losses.  This 
is  again  consistent  with  a  number  of  experimental  observations 
[4,48,50];  the  reference  exchange  current  densities  for  borohydride 
oxidation  (Table  1 )  and  hydrogen  evolution  ( ~  1CT4  A  m-2  in  NaOH 
at  room  temperature  [40])  on  Pt  are  much  higher  than  the  reference 
exchange  current  density  for  oxygen  reduction  on  Pt  (Table  1 ),  which 
partly  explains  these  observations.  It  should  be  pointed  out  that  the 
OCV  in  Fig.  3  is  between  0.2  and  0.1  V  higher  than  is  typically  found 
for  DBFC  employing  Pt/C  anodes  [6,7],  This  discrepancy  is  probably 
due  to  borohydride  crossover  to  the  cathode  (particularly  at  low 
current  densities)  and  the  establishment  of  a  mixed  cathode 
potential  [20,21],  which  is  not  incorporated  in  the  model. 


The  current  densities  associated  with  direct  borohydride 
oxidation  (4)  and  hydrogen  evolution/oxidation  (3)  are  shown  in 
Fig.  4,  together  with  the  mixed  anode  potential  and  mixed  anode 
open-circuit  potential.  The  borohydryide  current  density  (equal  in 
magnitude  to  the  hydrogen  evolution  current  density  at  the  OCV)  is 
non-zero  at  the  OCV  (Fig.  4(a)),  leading  to  the  establishment  of 
a  mixed  potential,  £mixed  of  around-0.87  V  (Fig.  4(b)).  This  potential 
is  considerably  more  positive  than  the  standard  potential  for 
reaction  (4),  as  is  invariably  the  case  in  practice.  At  an  applied 
current  density  of  around  125  mA  cm-2,  the  borohydride  and 
hydrogen  current  densities  are  equal  and  above  this  value  of  the 
applied  current  density  the  hydrogen  oxidation  current  density  is 
actually  higher  than  the  borohydride  oxidation  current  density.  It  is 
interesting  to  note  that  in  this  case  (i.e.,  for  the  given  kinetic 
parameters)  the  borohydride  current  density  remains  almost 
constant  as  the  applied  current  is  increased.  The  sum  of  the  two 
current  densities  (referred  to  the  same  geometric  area)  must  equal 
the  total  applied  current  density;  the  applied  and  hydrogen  current 
densities  differ  approximately  by  a  constant  value,  which  repre¬ 
sents  the  borohydride  current. 

The  mixed  open-circuit  anode  potential,  S^ixed 's  als°  shown  in 
Fig.  4(b).  This  value  represents  the  potential  at  which  the  net 
current  is  zero.  As  the  current  density  is  increased,  the  mixed  open- 
circuit  potential  increases,  reflecting  the  changes  in  the  BH4  surface 
concentration  and  the  adsorbed  species  coverages  at  steady  state. 
The  latter  are  shown  in  Fig.  5.  The  surface  coverage  of  hydrogen  is 
around  0.27  under  open-circuit  conditions  (Fig.  5(a)),  while  the 
coverage  of  BH4  is  considerably  lower.  This  result  was  robust  to 
order-of-magnitude  changes  in  the  value  of  the  reaction  constant 
k  bh»  for  the  given  values  of  the  hydrogen  evolution  kinetic 
constants  and  the  borohydride  exchange  current  density.  Outside 
this  range  of  values  for  the  adsorption  constant  k  BH,  the  results 
were  not  physically  meaningful  (given  the  highly  nonlinear  nature 
of  the  kinetic  model,  it  is  not  surprising  that  the  character  of  the 
solutions  is  sensitive  to  changes  in  the  kinetic  parameters).  The 
results  suggest,  therefore,  that  only  a  small  surface  coverage  of  BFff 
is  required  to  sustain  the  borohydride  oxidation  reaction,  given  the 
higher  overpotential  for  reaction  (4),  tjbh,  compared  to  that  for 
reaction  (3),  tjh.  At  a  zero  current  density,  the  BH4  surface 
concentration  (Fig.  5(b))  is  appreciably  lower  than  the  inlet  value, 
as  a  result  of  the  high  borohydride  current  density  (BH4 
consumption)  under  open-circuit  conditions  (Fig.  4(a)). 

Fig.  6  illustrates  the  relationships  between  the  applied  current 
density  and  (a)  the  cell  voltage  and  (b)  the  hydrogen  evolution/ 
oxidation  current  densities  for  different  BH4  concentrations.  The 


currents  and  (b)  the  mixed  anode  potentials  vs.  the  applied  (total)  current  density  for  the  base-case  parameters  in  Tables  1—4.  The  anode 


Fig.  4.  (a)  The  borohydride  and  hydrogen  ( 
catalyst  is  Pt/C. 
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Fig.  5.  (a)  The  borohydride  ion  and  hydrogen  adatom  surface  coverages  and  (b)  the  borohydride  ion  surface  concentration  vs.  the  applied  (total)  current  density  for  the  base-case 
parameters  in  Tables  1-4.  The  anode  catalyst  is  Pt/C. 


Fig.  6.  The  effect  of  varying  the  borohydride  ion  concentration  (in  2  M  NaOH)  on  (a)  the  cell  voltage  and  (b)  the  hydrogen  evolution/oxidation  current  density  (expressed  per  unit 
geometric  area  of  the  electrode)  vs.  the  applied  (total)  current  density.  The  other  parameter  values  are  given  in  Tables  1-4.  The  anode  catalyst  is  Pt/C. 


performance  improves  as  the  concentration  is  increased  (Fig.  6(a)), 
which  has  been  observed  consistently  in  experimental  cell  studies 
on  different  catalysts,  including  Pt  [21,32,33,51—53].  As  the  boro¬ 
hydride  ion  concentration  is  increased,  the  magnitude  of  the 
hydrogen  evolution  current  under  open-circuit  conditions 
increases  (Fig.  6(b)).  This  is  caused  by  an  increase  in  the  surface 
coverage  of  BH4  (by  a  factor  of  approximately  5  when  the  BH4 
concentration  is  increased  from  0.5  M  to  1.5  M),  while  the  change  in 
the  mixed  anode  potential  is  small  (approximately  0.02  V  more 
negative  at  the  OCV  between  0.5  M  and  1.5  M).  At  low  current 
densities,  the  performance  is  often  observed  to  deteriorate  as  the 
BH4  concentration  is  increased  [51].  Under  these  conditions,  the 
rate  of  BH4  crossover  increases,  which  lowers  the  (mixed)  potential 
in  the  cathode.  This  phenomenon  is  not,  however,  included  in  the 
model. 

Fig.  7  shows  the  cell  voltage  and  anode  polarization  for  different 
anolyte  flow  rates.  It  is  well  known  that  increasing  the  flow  rate 
(within  a  practical  range)  can  provide  relatively  small  improve¬ 
ments  in  performance  [11,15,21,52],  precisely  as  shown  in  Fig.  7.  At 
a  fixed  current  density,  a  faster  flow  rate  ensures  a  higher 
concentration  of  the  BH4  ions  in  the  anode,  which  lowers  the 
anode  polarization.  The  latter  result,  shown  in  Fig.  7,  agrees  with 
the  anode  polarization  measurements  taken  by  Duteanu  et  al.  [21  ]. 


polarization  vs.  the  applied  (total)  current  density.  The  other  parameter  values  are 
given  in  Tables  1— t.  The  anode  catalyst  is  Pt/C. 
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Fig.  8.  The  effect  of  varying  the  02  concentration  (pressure)  on  the  cell  voltage  and  the 
cathode  polarization  vs.  the  applied  (total)  current  density.  The  other  parameter  values 
are  given  in  Tables  l^t.  The  anode  and  cathode  catalysts  are  Pt/C. 


Varying  the  oxygen  concentration  in  the  cathode  inlet  stream 
can  have  a  greater  impact  on  DBFC  performance  than  varying  either 
the  flow  rate  or  the  BH4  concentration  (within  practical  ranges) 
[11,52],  Fig.  8  shows  the  cell  voltage  curves  and  combined  cathode 
activation/mass-transport  losses  for  cathode  air  pressures  of  1,  0.5 
and  0.25  bar  (assuming  21%  O2  by  volume).  There  is  a  significant 
improvement  in  performance  between  the  highest  and  lowest 
pressures;  the  peak  power  density  is  around  65%  higher  for  1  bar 
(the  peak  power  density  is  attained  at  95  mA  cm  2  for  1  bar  and  at 
69  mA  cm-2  for  0.25  bar).  The  same  degree  of  improvement  (ca. 
65%)  was  recorded  by  Celic  et  al.  [51]  when  air  was  replaced  with 
pure  O2,  i.e.,  a  factor  of  five  increase  in  the  concentration  of  O2.  This 
enhanced  performance  is  due  to  improved  mass  transport  to  the 
catalyst  sites  through  the  diffusion  boundary  layer  in  the  CCL 
(ionomer),  leading  to  a  higher  surface  concentration  of  O2.  The 
effect  of  the  higher  surface  concentration  on  the  cathode  over¬ 
potential  at  a  fixed  current  density,  which  is  the  main  factor  in  the 
performance  improvement,  is  evident  from  Fig.  8. 

It  is  worthwhile  examining  the  individual  components  of  the 
ohmic  losses,  the  most  significant  of  which  are  shown  in  Fig.  9(a) 
for  the  base-case  parameters  in  Tables  1—4.  The  losses  through  the 
electron  conducting  components  are  small  and  are  not,  therefore, 


Fig.  9.  (a)  The  main  components  of  the  total  ohmic  losses  vs.  the  applied  (total)  current  density  based  on  the  parameter  values  given  in  Tables  1  —4.  The  anode  catalyst  is  Pt/C.  Losses 
through  the  electron  conducting  components  were  small  and  are  not  included,  (b)  The  effect  on  the  voltage  drops  vs.  the  applied  (total)  current  density  of  increasing  the  CCL 
ionomer  volume  fraction  to  0.3,  with  all  other  parameters  fixed. 
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applied  (total)  current  density  in  the  case  of  a  Ni/C  anode  catalyst  The  parameter 
values  are  as  in  Tables  1^4  unless  otherwise  specified  in  the  text. 

included  in  this  figure.  For  the  given  diffusion  coefficients  and  inlet 
concentrations,  the  ionic  conductivity  of  the  anolyte  is  around 
10  S  m_1,  while  the  membrane/ionomer  conductivity  was  assumed 
to  take  a  value  of  0.5  S  m_1.  Therefore,  for  the  given  volume  frac¬ 
tions  in  Table  3,  the  effective  ACL  ionic  conductivity,  ff^cl,  defined 
in  equation  (32),  is  approximately  equal  to  the  effective  anolyte 
conductivity,  also  given  in  equation  (32).  The  ACL  losses  are  rela¬ 
tively  small  compared  to  the  equivalent  losses  in  the  CCL,  which  is 
a  consequence  of  the  parallel  pathways  for  ion  migration  in  the  ACL 
and  the  high  solution  conductivity  in  comparison  to  the  ionomer 
conductivity.  The  CCL  losses  are  strongly  influenced  by  the  volume 
of  ionomer,  and  can  be  reduced  significantly  by  increasing  this 
volume  fraction,  without  noticeably  increasing  the  mass-transport 
resistance  in  the  cathode,  as  demonstrated  in  Fig.  9(b). 

The  conductivity  of  the  membrane  is  an  important  parameter  of 
the  cell  and  can  be  altered,  at  a  fixed  temperature,  by  the  NaOFI  and 
water  concentrations  in  the  anolyte  (to  which  the  membrane  is 
exposed  on  the  anode  side)  and  the  relative  humidity  of  the  air  (to 
which  the  membrane  is  exposed  on  the  cathode  side)  [11  ].  The  type 
of  membrane/ionomer  also  influences  the  performance.  For 
example,  Ma  et  al.  demonstrated  a  substantial  uplift  in  performance 
by  replacing  a  Nation®  212  membrane/ionomer  pair  with  a  cross- 
linked  chitosan  membrane  (CCS)  in  combination  with  a  chitosan 


Fig.  12.  (a)  The  borohydride  and  hydrogen  currents  and  (b)  the  mixed  anode  potentials 
values  are  as  in  Tables  1-4  unless  otherwise  specified  in  the  text 


chemical  hydrogel  as  the  anode  binder;  the  conductivity  of  the  CCS 
equilibrated  with  10%  aqueous  NaOH  at  room  temperature  was 
recorded  11  S  m-1,  compared  to  0.74  S  nr  1  for  the  Nafion®  212 
membrane.  Fig.  10(a)  shows  the  predicted  ohmic  losses  when  the 
membrane  and  ionomer  conductivities  are  increased  from 
0.5  S  m-1  to  just  1  S  m-1,  with  all  other  parameters  as  in  Tables  1—4. 
Both  the  membrane  and  CCL  losses  are  halved,  the  latter  of  which 
leads  to  a  dramatic  reduction  in  the  total  ohmic  losses  (given  the 
small  volume  of  ionomer  in  the  CCL),  with  an  equal  reduction  in  the 
total  voltage  losses. 

Significant  variations  in  performance  have  been  observed 
between  cells  employing  different  Nafion®  membranes  [48,53],  The 
ohmic  loss  breakdown  in  Fig.  10(b)  shows  that  an  increase  in  the 
membrane  thickness  from  50  pm  to  175  pm  dramatically  increases 
the  membrane  losses  (the  peak  power  density  decreases  by  29%). 
The  thickness  values  of  50  pm  and  175  pm  correspond  to  the  dry 
thicknesses  of  Nafion®  112  (or  212)  and  Nafion®  117,  respectively, 
both  of  which  were  studied  by  Liu  et  al.  [48],  The  authors’  results 
showed  that  the  peak  power  density  decreased  by  14%  when 
Nafion®  112  was  replaced  with  Nafion®  117,  which  is  roughly  50%  of 
the  value  predicted  by  the  model.  If,  on  the  other  hand,  the 
membrane/ionomer  conductivity  is  increased  to  1  S  m-1,  the 
agreement  is  extremely  good  (both  around  14%  decrease  with 
Nafion®  117).  These  results  suggest  that  the  marked  difference  in 
performance  between  cells  employing  different  Nafion® 
membranes  is  largely  a  result  of  their  different  thicknesses,  rather 
than  any  major  differences  in  their  ion  transport  properties. 

3.2.  Performance  using  a  nickel  anode  catalyst 

To  simulate  the  performance  when  Ni  is  used  as  the  anode 
catalyst,  the  values  for  the  rate  constants  listed  in  Table  5  were 
used,  while  other  parameter  values  were  as  given  in  Tables  1  —4.  A 
value  of  the  exchange  current  density  for  the  borohydride  oxidation 
reaction  on  Ni,  jBH.ref  is  not  available  from  the  literature.  Cyclic 
voltammetry  results  [54,55]  suggest  that  the  exchange  current 
density  is  significantly  lower  than  that  on  Pt/C;  (/elikkan  et  al. 
concluded  that  the  activity  of  Ni  towards  borohydride  oxidation  is 
extremely  limited.  The  value  in  Table  1  was,  therefore,  divided  by  an 
arbitrary  factor  k,  which  was  varied  between  1  and  106.  It  was  found 
that  the  character  of  the  solutions  did  not  change  within  this  range; 
the  differences  were  seen  in  the  OCV  and  the  H  surface  coverage  at 
OCV,  and  will  be  discussed  below.  The  results  depicted  in 
Figs.  11—13  are  based  on  a  value  of  k  =  1000.  Likewise,  the  unknown 
adsorption  constant  for  reaction  (5),  k  BH  was  arbitrarily  chosen  as 


vs.  the  applied  (total)  current  density  in  the  case  of  a  Ni/C  anode  catalyst.  The  parameter 
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Fig.  13.  (a)  The  borohydride  ion  and  hydrogen  adatom  surface  coverages  and  (b)  the  borohydride  ion  and  hydrogen  molecule  surface  concentrations  vs.  the  applied  (total)  current 
density  in  the  case  of  a  Ni/C  anode  catalyst.  The  parameter  values  are  as  in  Tables  unless  otherwise  specified  in  the  text. 


one-tenth  of  the  value  for  Pt/C  in  Table  1.  This  parameter  was  also 
varied  as  discussed  below. 

Fig.  11  shows  the  cell  voltage  at  different  applied  current 
densities  and  the  individual  components  of  the  voltage  loss  for  the 
Ni  anode  cell.  The  equivalent  cell  voltage  curve  for  Pt/C  (taken  from 
Fig.  3)  is  included  for  comparison.  Comparing  the  cell  voltage 
curves  for  Ni  and  Pt/C  shows  that  on  Ni  the  OCV  is  around  0.13  V 
higher,  while  the  limiting  current  density  at  0  V  is  around 
45  mA  cm-2  lower.  There  are  few  direct  comparisons  between  Ni 
and  Pt  electrodes.  Wang  et  al.  [33]  found  a  difference  of  around 
0.13  V  in  the  OCV  between  Ni  powder  and  20  wt.%  Pt/C  anodes  (at 
296  I<  and  with  0.5  M  NaBFU),  with  a  higher  value  when  Ni  was 
used.  Liu  et  al.  [56]  recorded  a  difference  of  around  0.17  V  (at  303  K 
and  with  0.52  M  NaBFLj),  with  the  higher  value  again  for  Ni  powder 
vs.  5  wt.%  Pt/C.  Cheng  et  al.  [55],  on  the  other  hand,  found  that  the 
OCV  obtained  with  42  wt.%  carbon-supported  Ni  was  approxi¬ 
mately  0.03  V  lower  than  on  43  wt.%  Pt/C.  The  authors  also  found 
that  the  peak  power  density  with  Ni/C  was  21%  lower  than  for  Pt/C, 
which  compares  to  an  equivalent  decrease  of  17%  in  the  case 
depicted  in  Fig.  11.  Neither  Wang  et  al.  nor  Liu  et  al.  generated  full 
polarization  curves,  so  comparisons  cannot  be  made  to  their  data.  A 
discussion  on  the  aforementioned  trends  is  postponed  until  other 
aspects  of  the  simulation  results  are  discussed. 

Fig.  12  shows  the  borohydride  and  hydrogen  current  densities 
and  the  mixed  anode  potential  for  the  Ni  anode.  Comparing  this 
figure  with  Fig.  4  reveals  that  the  hydrogen  evolution  current 
density  at  the  OCV  is  substantially  lower  for  the  Ni  anode 
(Fig.  12(a))  and,  correspondingly,  the  mixed  anode  potential  is  more 
negative  (Fig.  12(b)),  which  leads  to  the  higher  OCV.  As  is  evident 
from  the  rate  constants  in  Tables  1  and  5,  Ni  is  less  active  than  Pt/C 
towards  hydrogen  evolution  and  borohydride  oxidation.  Both  the 
borohydride  and  hydrogen  current  densities  at  the  OCV  are, 
therefore,  lower  than  for  Pt/C.  The  surface  coverages  of  both  H  and 
BHJ  at  the  OCV,  shown  in  Fig.  13(a),  are  higher  than  the  coverages 
on  Pt/C  in  Fig.  4(a).  These  increases  result  from  the  slower  boro¬ 
hydride  and  hydrogen  reactions  on  Ni. 

The  direction  in  which  the  mixed  anode  potential  at  the  OCV 
moves  (towards  £{]H  or  E^)  as  the  reaction  constants  are  varied 
depends  on  the  relative  sizes  of  the  exchange  current  densities  for 
the  borohydride  and  hydrogen  reactions.  As  jBH.ref  is  decreased 
(with  all  other  parameters  fixed)  by  increasing  the  fitting  param¬ 
eter  k,  the  mixed  potential  at  a  fixed  current  density  moves  closer  to 
the  standard  potential  for  hydrogen  evolution/oxidation.  The 
borohydride  overpotential  t/bh  increases  in  order  to  maintain 


a  balance  between  the  hydrogen  and  borohydride  current  densities 
(see  equations  (16)  and  (18)),  which  leads  to  a  smaller  hydrogen 
overpotential  tjh  (the  mixed  anode  potential  is  less  negative)  and 
concomitantly  a  smaller  hydrogen  evolution  current  density.  The 
lower  cell  voltage  values  lead  to  a  lower  peak  power  density.  For 
k  =  2  x  104,  the  OCV  is  approximately  0.02  V  higher  and  the  peak 
power  density  is  approximately  21%  lower  than  on  Pt/C,  which  is 
very  close  to  the  result  of  Cheng  et  al.  for  a  cell  employing  Ni/C  [55], 
The  differences  between  the  results  of  Wang  et  al.  and  Liu  et  al. 
[33,56]  and  those  of  Cheng  et  al.  [55]  could,  therefore,  be  explained 
by  the  different  exchange  current  densities  for  the  Ni/C  and  the  Ni- 
powder  catalysts  used. 

The  adsorption  constant  for  reaction  (5),  kBH  was  arbitrarily 
chosen  for  the  Ni  simulation.  Decreasing  the  value  of  kBH  was 
found  to  decrease  the  cell  voltage  at  a  fixed  current  density  by 
virtue  of  a  lower  surface  coverage  of  BH4 ,  requiring  a  larger  over¬ 
potential  for  borohydride  oxidation,  ijbh  in  order  to  maintain  charge 
conservation  (equal  borohydride  and  hydrogen  current  densities). 
The  resultant  decrease  in  the  hydrogen  overpotential,  tjh  decreases 
the  hydrogen  evolution  current  density.  The  basic  character  of  the 
solutions  (as  described  above)  is  not  altered  unless  values  of  k  Bh 
smaller  than  4  x  10-7  mol  m2  s-1  are  used,  in  which  case  the 
borohydride  surface  coverage  is  limiting.  The  surface  concentration 
of  BH4  reaches  a  limiting  value,  along  with  the  borohydride  and 
hydrogen  overpotentials,  leading  to  a  characteristic  bend  in  the 
polarization  curve.  There  is  no  concrete  example  of  such  behaviour 
in  the  literature,  so  these  low  values  of  k  Bh  were  not  considered  to 
be  physically  meaningful. 

4.  Conclusions 

In  this  paper,  a  framework  for  the  detailed  modelling  of  direct 
borohydride  fuel  cells  has  been  developed.  Comparisons  to 
experimental  data  from  the  open  literature  have  demonstrated 
that  the  model  is  able  to  capture  the  correct  qualitative  trends 
with  respect  to  variations  in  several  key  parameters,  including 
the  fuel  and  oxidant  concentrations  and  the  membrane  proper¬ 
ties.  A  comparison  between  the  simulated  performance  for  Pt 
and  Ni  anode  catalyst  cells  has  revealed  a  number  of  key  differ¬ 
ences.  These  differences  were  explained  in  the  context  of  the 
trends  that  have  been  observed  in  experimental  studies.  It 
appears  that  the  hydrogen  evolution  (water  reduction)  current 
on  Ni  is  considerably  lower  at  the  OCV  (given  the  lower  activity 
of  Ni  towards  this  reaction),  with  a  correspondingly  more 
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negative  mixed  anode  potential,  which  ultimately  leads  to 
a  higher  OCV.  The  results  also  suggest  that  the  OCV  decreases  as 
the  exchange  current  density  for  borohydride  oxidation 
decreases,  which  may  explain  the  differences  between 
Ni-powder  and  Ni/C  catalysts. 

The  anode  mechanism  is  highly  complex;  it  depends  on  the 
concentration  of  sodium  borohydride,  the  concentration  of  the 
alkali  salt  (NaOH  or  KOH  typically),  the  temperature  and  the  cata¬ 
lyst.  The  anode  kinetic  mechanism  adopted  in  this  paper  incorpo¬ 
rates  borohydride  oxidation,  hydrogen  evolution  through 
hydrolysis,  hydrogen  evolution  through  water  reduction  and 
hydrogen  oxidation.  The  role  of  intermediates  such  as  BFIs/OH)^ 
and  the  adsorption  of  OH-  may,  on  the  other  hand,  play  important 
roles.  It  is  important,  therefore,  to  develop  this  aspect  of  the  model 
further,  which  requires  values  for  the  rate  constants  pertaining  to 
the  various  reaction  mechanisms  proposed  in  the  literature  (for 
each  catalyst  materials  of  interest). 

Other  important  parameters  that  have  not  been  measured 
systematically  relate  to  the  transport  of  species  in  the  anolyte 
and  the  transport  of  Na+  in  the  membrane.  The  simulations 
confirm  that  ohmic  losses  represent  a  considerable  portion  of  the 
total  cell  losses.  The  Na+  conductivities  of  the  perfluorosulfonic- 
acid  and  other  ion-exchange  membranes  used  in  DBFCs  depend 
on  the  operating  environment.  The  conductivity  values  can  vary 
significantly  between  membrane  types  and  between  operating 
conditions.  In  order  to  fully  exploit  models  for  design  and  char¬ 
acterisation,  it  is  particularly  important  to  develop  relationships 
between  the  Na+  conductivity  and  the  membrane  water  content/ 
temperature  of  the  environment. 

The  aim  of  this  work  was  to  develop  a  first  detailed  model  of 
a  DBFC.  Aspects  for  consideration  in  extending  the  present  model 
include  spatial  variations  in  the  reactant  concentrations  and 
potentials,  heat  management  and  dynamic  performance.  The 
framework  developed  allows  these  aspects  of  performance  to  be 
incorporated  in  a  relatively  straightforward  manner. 


Appendix 


Charge  conservation  and  electrical  neutrality  in  the  catalyst 
layers  demand  that  the  electronic  and  ionic  current  densities,  js  and 
je,  satisfy  (neglecting  double  layer  effects): 


Integrating  this  equation  and  applying  the  boundary  condition  at 
either  the  membrane/catalyst-layer  interface  (zero  electronic 
current  density)  or  the  diffusion-layer/catalyst-layer  interface  (zero 
ionic  current  density)  yields: 

je  +  js  =  ./applied  (A2) 

The  total  (applied)  current  density  (referred  to  the  geometric 
surface  area  of  the  electrode)  is  shared  between  the  electronic  and 
ionic  conducting  phases.  The  individual  components  je  and  js  are 
given  in  terms  of  the  potentials  0e  and  0S  by  ohm’s  law,  as  in 
equations  (33)— (36).  They  can  also  be  related  to  the  transfer 
(Faradaic)  current  densities,  jf  in  the  electrodes  by  the  following 
charge  balances  (consistent  with  equation  (Al))  [22]: 


(A3) 


where  as  and  <re  are  the  relevant  effective  electronic  and  ionic 
conductivities  and  s  is  the  specific  surface  area  of  the  electrode 
(active  surface  area  per  unit  volume).  The  transfer  current  density  is 
given  byjo2  in  the  CCL  and  by j'bh  +Jv  +Jh  in  the  ACL  (equations  (15) 


and  (17)).  From  their  definitions,  the  specific  surface  area  can  be 
related  to  the  roughness  factor,  r,  and  the  thickness  of  the  catalyst 
layer,  /,  as  s  =  r//.  Therefore: 


d2</>e 
Te  dx2 


d  0S  -/'applied 

s  dx2^ - 


(A4) 


Integrating  the  second  of  these  equations  twice  from  the  liquid 
diffusion  layer/ACL  interface,  x  =  Xi,  to  some  arbitrary  x  in  the  ACL, 
applying  the  boundary  condition  (7s(d0s/dx)  =  ./applied  at  x  —  xi, 


Js  =  HvT  =  - 


(A5) 


The  average  value  of  js  is  found  by  integrating  the  right  hand  side  of 
equation  (A5)  w.r.t.  x  between  x  =  xi  and  the  membrane/ACL 
interface,  x  =  X2,  which  gives  a  value  of  jappiied/2,  as  used  in  the 
formulae  in  Section  2.4.  Integrating  equation  (A5)  across  the  ACL 
yields  the  equivalent  formula: 

A0  =  0(x2)-0(xt)  =  (A6) 


This  analysis  can  be  extended  to  the  ionic  potential  in  the  ACL  and 
to  the  electronic  and  ionic  potentials  in  the  CCL. 
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